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La presente invention a trait a une methode pour mesurer les similarites locales" 
entre plusieurs cubes de traces sismiques 3D avant sommation obtenus d'un volume d'une 
zone souterraine, ou a la suite de campagnes de prospection repetitives (4D), par exemple. 
II s'agit done d'une raesure locale de coherence qui traduit, en premier lieu, la 
ressemblance d'un cube sismique a 1'autre, tout en prenant en compte la ressemblance 
locale au sein d'un meme cube. 

ETAT DE LA TECHNIQUE 

La notion de coherence proprement dite est une notion relativement' recente. La 
problematique, jusqu'a present, a ete de mettre au point un outil revelant les changements 
stratigraphiques ou structuraux (notamment les failles) a partir des mesures sismiques, et^ 
d'obtenir ainsi une information volumique sur ces changements. Le fondement de toutes 
les methodes developpees depuis moins de dix ans consiste a definir une dissimilarity 
locale de trace a trace. / v s 

Un premier algorithms decrit par: Bahorich, M., and Farmer, S. (1995), «The 
coherence cube », The Leading Edge, 14, 10, 1053-1058, consiste a calculer la correlation > 
croisee entre chaque trace d'un cube sismique avec ses deux voisines en ligne, avec ses 
deux voisines en CDP (couverture multiple), puis a combiner les deux resultats, apres les 
avoir normalises par 1'energie des traces. La coherence n'est estimee qu'a partir de trois 
traces, ce qui rend le calcul tres rapide mais aussi peu robuste si les donnees sont bruitees. 

Suivant un autre algorithme decrit par Marfurt, K. J. 5 et al. (1998), « 3-D Seismic 
Attributes Using a Semblance-based Coherency Algorithm », Geophysics, 63, 1150-1165, 
le calcul de la coherence est base sur un calcul de semblance locale mettant en jeu 
davantage de traces, ce qui rend le resultat plus robuste au bruit. 
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Suivant un autre algorithme decrit par Gersztenkorn, A., et Marfurt, K. J. (1999), 
« Eigenstructure based Coherence Computations as an Aid to 3-D Structural and 
Stratigraphic Mapping ». Geophysics, 64, 1468-1479, le calcul de la coherence est base sur 
une decomposition en valeurs propres : une fenetre d'analyse definie en lignes, CDP et 
5 temps est extraite du cube sismique, la matrice de covariance des traces sismiques est 
constitute, et la plus grande valeur propre de cette matrice est calculee. La valeur de 
coherence correspond alors au rapport entre cette valeur propre et la somme de toutes les 
valeurs propres de la matrice de covariance, ou trace de la matrice de covariance, qui n'est 
autre que la variance totale des traces sismiques de la fenetre d'analyse. 

.10 Toutes ces approches presentent certaines limites. En particulier, une limitation 

majeure est qu'elles ne sont pas applicables a r analyse de donnees sismiques multicubes. 

En effet, le but de ces differents attributs de coherence est plutot de cartographier 
les anomalies stratigraphiques : ils ne permettent done pas d'evaluer la coherence soit 
calendaire (4D) soit AVO (pour "Amplitude Versus Offset"). II n'existe, a ce jour, a notre 
15 connaissance, aucun algorithme permettant de determiner de tels attributs. 

L' Analyse en Composantes Principales Generalised (ci-apres ACPG) est un outil 
connu permettant de mettre en evidence une eventuelle redondance d'information entre des 
groupes d' attributs sismiques : elle peut etre adaptee pour definir une mesure de 
ressemblance locale des donnees sismiques, et ce d'un cube a T autre, en analysant un 
20 voisinage autour d'un point courant, la notion de groupe d' attributs etant reliee aux 
differentes campagnes dans le temps ou aux differentes campagnes sismiques avant 
sommation, par exemple. 

Cette technique a ete mise en ceuvre dans la methode decrite dans la demande de 
brevet FR 02/1 1 .200 du demandeur, pour compacter et filtrer des evenements sismiques lus 

25 sur des traces sismiques "multicubes", avec repartition de ces evenements en families 
correspondant chacune a une signification physique particuliere : cube de donnees iso 
deport ou iso angle d'incidence, cubes de parametres elastiques issus d'une inversion 
strati graphique jointe, etc. dans le but d'en extraire des informations sur la nature du sous- 
sol. Cette methode comporte essentiellement la formation par combinaison des variables 

30 sismiques, de variables synthetiques en nombre tres inferieur, que Ton obtient par 
construction d'une base vectorielle orthogonale dans chacun des ensembles d'analyse 
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constitues par les donnees de chacune des families, d'ou decoule la formation d'une base 
vectorielle orthonormee decrivant ces ensembles d' analyse, et l'utilisation de cette base 
vectorielle orthonormee (nouveaux attributs) pour filtrer et decrire les dits evenements 
sismiques. 

5 La methode selon l'invention 

La methode selon l'invention propose une mesure de la similarity locale entre 
plusieurs cubes de donnees sismiques 3D avant sommation ou sismiques 4D (repetees dans 
]e temps). La methode comporte les etapes suivantes : 

a) en chaque point du volume etudie et caracterise par plusieurs cubes sismiques, on 
10 extrait un voisinage volumique, centre sur ce point (point courant) et compose d'un 

ensemble de traces sismiques en nombre limite ; ainsi, chaque point courant est caracterise 
par autant de groupes d' attributs sismiques qu'il y a de cubes analyses ; 

b) on applique la technique d'analyse ACPG a ces groupes d' attributs sismiques 
extraits de chacun des cubes sismiques dans le voisinage volumique du point courant pour 

15 former des variables synthetiques ; 

c) on determine une valeur de coherence a partir des variables synthetiques 
extraites, que Ton affecte au point courant ; 

d) on reitere les etapes a) a c) pour chaque point ; et 

e) on regroupe Tensemble des valeurs de coherence en un cube de coherence. 

20 Les valeurs comprises dans ce cube de coherence donnent le degre de ressemblance 

locale recherche entre les cubes de donnees sismiques. 

Les projections des variables synthetiques sur les differents cubes dans le voisinage 
du point courant represented une part de Finformation du groupe correspondant. Cette part 
d' information ou de variance est connue. Des lors> plusieurs approches peuvent etre 
25 envisagees pour le calcul de l'attribut de coherence a partir des valeurs de correlation 
calculees entre les variables synthetiques et leurs projections sur les cubes dans le 
voisinage du point courant. 
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Suivant un mode de mise en oeuvre, pour chaque point, on prend comme valeur de 
coherence, la valeur moyenne des carres des correlations entre les variables synthetiques et 
leurs projections sur les cubes dans le voisinage du point courant, sur un nornbre limite k 
des dites variables synthetiques. 

5 La valeur de k est determinee, par exemple, comme etant le plus petit nombre de 

variables synthetiques permettant d'atteindre un seuil de variance explique par les 
projections des variables synthetiques sur chacun des cubes, ce seuil etant prealablement 
choisi. 

Suivant un autre mode de mise en oeuvre, on choisit un nombre de variables 
10 synthetiques selon leurs correlations avec les groupes d'attributs associes au voisinage 
volumique du point courant. La valeur de coherence attribuee au point courant est egale a 
la somme ponderee des carres des correlations entre les variables synthetiques considerees 
et leurs projections sur les cubes dans le voisinage du point courant. 

Pour une valeur de correlation, on choisit par exemple comme valeur de 
15 ponderation, le pourcentage de variance explique par la projection de la variable 
synthetique sur le groupe correspondant divise par la somme des variances de toutes les 
projections des variables synthetiques considerees sur le meme groupe. 

Suivant un autre mode de mise en ceuvre, on fixe un seuil sur le pourcentage de 
variance explique par les projections des variables synthetiques sur les cubes, dans le 
20 voisinage du point courant, qu'il faut prendre en compte. La valeur de coherence est alors 
egale a la somme ponderee des carres des correlations entre les variables synthetiques et 
leurs projections sur les cubes dans le voisinage du point courant, de sorte que le nombre 
de variables synthetiques pris en compte permette d'atteindre ce seuil. 

Pour une valeur de correlation definie, on choisit par exemple une valeur de 
25 ponderation egale a p (nombre de cubes) fois le seuil de variance fixe. 

Selon les cas, le voisinage volumique peut etre extrait de cubes de traces sismiques 
obtenus soit a Tissue d'une campagne sismique 3D, chacun correspondant a un meme angle 
d'incidence ou a un meme deport ou offset, soit a Tissue .de campagnes d'exploration 
sismiques successives de la zone. 
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Le voisinage volumique peut etre aussi extrait de cubes de residus obtenus soil a 
Tissue d'une inversion strati graphi que avant sommation soit de cubes de residus obtenus a 
Tissue d'une inversion strati graphique apres sommation.. L peut etre aussi extraits des 
cubes inverses (avant ou apres sommation) et des cubes de residus. 

La methode est particulierement avantageuse en ce qu'elle permet de definir un 
nouvel attribut mesurant une ressemblance locale entre plusieurs cubes sismiques extraits 
d'un voisinage autour d'un point. Elle permet de prendre en compte T aspect multicube des 
donnees sismiques et mesure davantage la variability d'un cube sismique a. F autre que la 
'variability au sein d'un meme cube. 

Presentation succingte des figures 

la figure 1 represente Textraction des cubes sismiques pour Tanalyse de coherence, 
dans le voisinage d'un point courant ; 

la figure 2 montre les projections de la variable synthetique Z® sur les groupes bet k ; 

.V 

la figure 3 represente les cubes sismiques (a) (b) et (c) obtenus a la suite de trois 
campagnes sismiques repetees dans le temps et le cube de coherence associe (premier 
mode de mise en oeuvre ou approche) - Fenetre temporelle en dehors du reservoir ; 

la figure 4 montre des lignes extraites du cube de coherence - (a) ligne 10, (b) ligne 20, 
(c) ligne 30, (d) ligne 40 ; 

la figure 5 montre un plan situe a 28ms sous le toit extrait des cubes des memes trois 
campagnes et du meme cube de coherence ; 

la figure 6 montre la ligne 10 extraite du cube de coherence calcule selon le premier 
mode de mise en oeuvre avec un seuil a 99% (a), selon le troisieme mode avec un seuil 
a 99% (b), a 90% (c), selon le.deuxieme mode de mise en oeuvre avec la premiere 
variable synthetique (d), les deux premieres variables synthetiques (e) ; 

la figure 7 montre des exemples de distributions des differences d' amplitude entre deux 
campagnes sismiques ; 

la figure 8 montre. les cubes sismiques associes aux trois campagnes successives et le 
cube de coherence associe, dans une fenetre temporelle au niveau du reservoir ; 
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la figure 9 montre les plans temporels extraits du cube de coherence calcule sur le 
reservoir ; 

- la figure 10 montre le plan temporel situe 12ms sous le toit d'un reservoir et Pattribut 
de coherence calcule avec la premiere variable synthetique (a), les deux premieres 

5 variables synthetiques (b), avec un seuil de variance a 90% (c),.a 95% (d), a 99% (e) ; 

- la figure 1 1 montre une vue 3D du cube de coherence obtenue avec les deux premieres 
variables synthetiques (seconde approche) - valeurs de coherence inferieures 
strictement a 0,8; 

- la figure 12 montre les cubes sismiques iso-angles 0°-6°, 12°-18°, 24°-30° et le cube de 
10 coherence associe ; 

- la figure 13 montre trois plans temporels situes (a) 4ms, (b) 10ms, (c) 16ms, sous le toit 
du reservoir et extraits du cube de coherence ; et 

- la figure 14 montre une ligne passant par un puits W2 extraite des cubes sismiques 0°- 
6°, 12°-18°, 24°-30° et du cube de coherence. 

15 Description detaillee 

La notion de coherence jusqu'ici a ete surtout appliquee a la recherche d'anomalies 
strati graphiques et les valeurs de coherence calculees a partir d'un seul cube de donnees 
sismiques, habituellement le cube apres sommation ou "stack". 

Avec la methode que Ton va decrire ci-apres, on forme un cube de coherence a 
20 partir de plusieurs cubes de donnees sismiques 3D (AVO ou 4D), indiquant en tout point le 
degre de ressemblance ou de dissimilarity locale des donnees sismiques, cube a cube, sur 
un voisinage volumique autour d'un point courant, et permettant done d' imager ce qui 
change ou ce qui ne change pas d'un cube a P autre. 

Comme on l'a rappele, TACPG est une technique permettant de mettre en evidence 
25 ce qui est commun et ce qui est different entre p groupes de variables ou d'attributs 
sismiques, et de determiner rapidement si tous les groupes sont lineairement identiques. 
Dans le cas du calcul d'un cube de coherence, il s'agit de realiser une mesure locale de la 
ressemblance (ou dissemblance) d'un cube sismique a 1' autre, tout en prenant aussi en 
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compte la ressemblance locale autour du point courant au sein d'un meme cube. 

Considerons p cubes de traces sismiques. Ces cubes de traces peuvent 
corresponds, par exemple, a des campagnes sismiques apres sommation, repetees dans le 
temps sur une meme zone geographique (cubes sismiques 4D), ou bien a des cubes 
5 sismiques 3D avant sommation iso angle ou iso-"offsets M (deports). 

Un voisinage volumique, centre sur une coordonnee (Ligne ; CDP (couverture 
multiple) ; temps ou profondeur) et compose d'un nombre limite de traces, est extrait de 
chacun des p cubes sismiques (figure 1). Le nombre de traces composant ce voisinage sera 
discute plus loin. On dispose ainsi de p ensembles de traces de meme dimension centres 
10 sur un point de memes coordonnees geographiques, et correspondant aux p cubes 
sismiques initiaux. 

Une ACPG est realisee sur les p ensembles ainsi extraits. Chaque ensemble ainsi 
extrait dans le voisinage du point courant correspond a un groupe d'attributs sismiques 
initiaux, ces attributs etant simplement, par exemple, les series des valeurs d'amglitudes 
15 correspondant aux differentes valeurs de la trace dans la fenetre temporelle etudiee. Le 
nombre total d'attributs est done egal a p fois la dimension verticale du voisinage 
considere. 

On calcule alors le carre de la correlation entre la variable synthetique Z 0) et sa 
projection sur un groupe d'attributs (figure 2). Le carre de la correlation correspond, en 
20 fait, au cosinus carre de Tangle 8 entre les deux vecteurs representant respectivement la 
variable synthetique et sa projection. II donne une indication sur le degre de proximite 
entre ces deux vecteurs, et done entre la variable synthetique Z® et le groupe 
correspondant : une valeur de 1 indiquera que la variable synthetique et sa projection sont 
confondues, alors qu'une valeur eloignee de 1 donnera une indication sur leur eloignement. 

25 Ainsi, lorsque tous les groupes d'attributs se ressemblent, les cosinus carres des 

angles entre tous les Z (i) et leurs projections sont egaux a 1. Dans le cas contraire ou la 
ressemblance est faible, les carres des correlations sont plus ou moins eloignes de 1 pour 
un certain nombre de Z (i) ; ils sont d'autant plus eloignes, pour nombre d'entre eux, que les 
groupes d'attributs seront differents. 
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Or, les projections de chacun des Z (i) sur les differents groupes representent une 
part de reformation du groupe correspondant. Cette part de variance peut etre connue et 
calculee. Des lors, plusieurs approches peuvent etre envisagees pour le calcul de I'attribut 
de coherence a parti r de ces valeurs de correlation. 

5 Premiere approche 

Une premiere approche, simple, consiste a calculer la valeur moyenne des carres 
des correlations sur un nombre k de Z® (k < p). Le nombre k est choisi de la fa?on 
suivante : 

(i) un seuil S sur la variance cumulee est fixe, par exemple 90% ; 

10 (ii) k est alors determine comme etant le plus petit nombre de variables synthetiques 

Z (j) permettant d'atteindre ce seuil. 

Dans ce cas, le nombre de variables synthetiques considerees dans le calcul des 
correlations est identique pour chaque groupe et le poids attribue a chaque correlation est le 
meme. 

15 c=^££p 2 (Z<»,ZP) 
Seconde approche 

Une seconde approche consiste a choisir le nombre de variables synthetiques 
selon leur correlation avec les groupes : en general, les premieres variables suffisent 
puisque, par principe, elles representent la partie de P information commune aux groupes. 

20 Une fois ce nombre fixe, et a la difference de la premiere approche, la somme, 

ponderee par les variances, des carres des correlations entre les consideres et leurs 
projections sur les groupes est calculee. Les carres des correlations entre un vecteur Z (j) et 
ses projections sur les differents groupes peuvent en effet Stre tous egaux a 1, alors que la 
part de variance expliquee est faible. 

25 Ponderer par la variance permet alors de prendre en compte la capacite de 

compaction des variables synthetiques extraites de 1'ACPG dans le calcul de la coherence, 
et d'eviter d'attribuer une valeur trop forte si, en realite, les cubes de traces etudies ne se 
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ressemblent que pour une faible part. Dans ce cas, le poids pij affecte a chaque correlation 
est egal a la variance expliquee par la projection de la variable synthetique~Z 0) sur le 
groupe i correspondant, divisee par la somme de toutes les variances. Cette 
« normalisation » assure que la somme des poids est bien egale a 1. 

i=l j=l 

Outre la difference de ponderation avec la premiere approche, il est & noter que la 
part de variance prise en compte dans chacun des groupes peut etre differente. 

Troisieme approche 

Enfin, une troisieme approche consiste, comme dans la premiere methode, a fixer 
10 un seuil sur la part de variance totale expliquee a prendre en compte. Mais cette fois-ci, 
pour chaque groupe i, le nombre ki de variables synthetiques Z (i) considerees sera 
strictement.le nombre permettant d'atteindre le seuil. Ainsi, ce nombre peut differer d'un 

K 

groupe a l'autre. La correlation « moyenne » sera estimee avec Fensemble des correlations 
elementaires des variables synthetiques necessaires pour chaque groupe. 

i=i j=i . ■ * 

Le poids pi j donne a chaque correlation est alors egal a la variance expliquee par la 
projection de la variable synthetique Z (i) sur le groupe i divise par p fois le seuil de 
variance choisi. Cette « normalisation » permet ainsi d'avoir une somme des poids egale a 

20 Deux parametres caracterisant la taille du voisinage d'analyse autour du point 

courant sont a determiner : le nombre de traces du voisinage et la dimension verticale (en 
temps ou en profondeur) des traces. Si un petit nombre de traces est pris en compte, par 
exemple neuf traces par voisinage, le resultat apparaitra spatialement plus bruite que si 
chaque voisinage est constitue de davantage de traces, par exemple 25. D'autre part, plus la 

25 dimension verticale sera importante, plus on peut s'attendre a ce que le resultat en 
coherence soit lisse verticalement. De plus, la variabilite pouvant augmenter, le seuil de 
variance a fixer dans le calcul de Tattribut de coherence selon la troisieme methode sera 
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different en fonction de la dimension verticale de la fenetre d' analyse. De meme, on peut 
s'attendre a ce que la capacit6.de compaction des variables synthetiques soit d'autant plus 
forte que la dimension de la fenetre est petite. 

EXEMPLES D' APPLICATION 

1 - Application a des donnees sismiques 4D 

La sismique repetee consiste a realiser plusieurs campagnes sismiques sur une 
meme zone geographique afin d'analyser et d' imager les changements pouvant intervenir 
au sein d'un reservoir apres sa mise en production. Calculer un attribut de coherence sur 
des donnees 4D a deux objectifs : 

1. indiquer plus precisement la reproductibilite du signal sismique en dehors 
du reservoir et controler ainsi le processus d'homogeneisation des amplitudes sismiques; 

2. indiquer ou et dans quelle mesure la reponse sismique varie au sein du 
reservoir et done aider a V interpretation de ces changements. 

On a utilise les traces sismiques de trois cubes apres sommation correspondant a 
trois campagnes sismiques 3D, desquels ont ete extraits trois cubes de 60 ras d'epaisseur 
chacun, situes environ 70ms au-dessus du reservoir, et trois cubes de 20 ms d'epaisseur 
situes au niveau du reservoir, 

L/ analyse des cubes en dehors du reservoir a pour but d'etudier dans quelle mesure 
le signal sismique se repete d'une campagne a 1'autre, alors que l'analyse des cubes 
sismiques situes au niveau du reservoir permet d'etudier les variations de la sismique au 
cours du temps induites par Fexploitation du reservoir. 

1-1 En dehors du reservoir 

Un attribut de coherence a tout d'abord ete calcule selon la premiere methode de 
calcul sur une partie situee largement au-dessus du reservoir (70ms) afin que les 
enregistrements sismiques ne soient pas influences par I'exploitation du reservoir. Le seuil 
de variance a ete fixe a 99%, permettant ainsi de prendre en compte la quasi-totalite de 
rinfomiation expliquee par les variables synthetiques extraites de 1'ACPG, mais aussi de 
ne pas prendre en compte des variables synthetiques expliquant une trop faible part de la 
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variance. La taille du voisinage du point courant utilise pour ]e calcul de la coherence est 
de 25 traces (un cube de 5 traces de cote centre au point courant) de 4ms chacune. La 
figure 3 montre les trois cubes sismiques correspondant aux trois campagnes, ainsi que le 
cube de coherence associe. Contrairement a ce qui pouvait etre attendu, ce dernier montre 
5 que les trois campagnes ne sont pas parfaitement coherentes puisque Ton trouve des 
valeurs inferieures a 0,7. 

Les trois campagnes sismiques semblent relativement bien coherentes sur les 
premieres 22 ms avec une majorite de valeurs superieures a 0,8 (figure 4). Au-dela, il 
existe, localement, davantage de zones presentant une faible valeur de coherence, avec une 
10 majorite de valeurs comprises entre 0,7 et 0,8 et, localement, des valeurs inferieures a 0,7. 

Ceci est illustre par la figure 5 qui montre le plan temporel, situe a +30ms en 
dessous du toit du cube, pour les trois campagnes sismiques ainsi que le cube de 
coherence. Les valeurs inferieures a 0,8 sont majoritaires et se repartissent sur Tensemble 
du plan temporel. Les sections sismiques des trois campagnes confirment bien ce manque 
15 de coherence 4D. 

Les cubes de coherence selon les deux autres methodes ont ete egalement calcules a 
partir des memes cubes sismiques. 

La figure 6 montre la ligne 10 extraite des cubes de coherences calcules selon la 
premiere methode avec un seuil fixe a 99% (a), selon la troisieme methode avec seuil a 
20 99% (b), a 90% (c), selon la deuxieme methode avec la premiere variable synthetique (d), 
les deux premieres variables synthetiques (e). 

Toutes les sections obtenues sont globalement assez proches. La section (c) montre 
des valeurs de coherence plus fortes que la section (b): la part de variance supplemental 
prise en compte semble done correspondre a une part d'information locale moins 
25 commune, tirant ainsi vers le bas la coherence. 

Les valeurs de coherence apparaissent un peu plus fortes lorsqu'elles sont 
ponderees par la variance, que lorsqu'une simple moyenne est calculee. La section (e) est 
similaire a la section (b) et la section (d) est similaire a la section (c) : il semble done que, 
dans la majorite des cas, localement, deux variables synthetiques suffisent pour resumer la 
30 totalite de V information. 
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La section (e) presente un peu plus de valeurs de faibles coherences que la section 
(d). De meme, les zones de tres fortes coherences (valeurs superieure a 0,9) sont un peu 
moins etendues dans le second cas. En revanche, quelques rares zones voient leur 
coherence legerement augmenter. Global ement, les resultats obtenus ne sont pas 

5 fondamentalement differents, bien qu'en ajoutant la seconde variable synthetique dans le 
calcul de r attribut de coherence, davantage de variance soit prise en compte. L'ajout de la 
seconde variable synthetique confirme done les ressemblances ou les dissemblances deja 
constatees avec une seule variable synthetique attribut. En conclusion, pour cette analyse 
menee hors reservoir, une seule variable synthetique peut suffire pour calculer F attribut de 

10 coherence. 

Les resultats ne sont pas detailles ici, mais il a ete verifie qu'en diminuant le 
nombre de traces definissant le voisinage (9 au lieu de 25), les cubes de coherence obtenus 
ont un aspect spatialement plus bruite. De meme, il a ete verifie qu'en augmentant la 
dimension verticale des traces sismiques, le cube de coherence obtenu etait lisse 
15 verticalement : dans ce cas, les tres faibles valeurs de coherence constatees sur la figure. 5 
sont legerement plus fortes. En prenant en compte deux ou trois variables synthetiques, ou 
en fixant un seuil de variance a 99%, il apparait moins de zones presentant de faibles 
valeurs de coherence. 

Quelle que soit la methode, il apparait que les cubes situes en dehors du reservoir 
20. ne sont pas totalement coherents : ceci peut etre du soit a un processus d'homogeneisation 
des amplitudes imparfait, soit a une certaine influence de l'exploitation du resei"voir sur les 
amplitudes. 

La figure 7 montre les distributions des differences d' amplitude entre deux 
campagnes successives a plusieurs annees d'intervalle, au sein de la fenetre temporelle 

25 etudiee. En cas de parfaite reproductibilite du signal, les valeurs mediane ou moyenne 
devraient etre centrees sur 0, et les distributions devraient etre peu etalees. Or, il apparait 
nettement que cette hypothese n'est exacte qu'entre 8 et 24 ms dans l'exemple considere. 
Ailleurs, les distributions oscillent autour de 0, avec une valeur mediane maximum atteinte 
aux alentours de 30 ms. Cette mesure globale des differences d'amplitude confirme bien le 

30 resultat plus local obtenu avec 1' attribut de coherence. 
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1-2 Dans le reservoir 

Un attribut de coherence a ensuite ete calcule au sein meme du reservoir selon la 
premiere methode. Le seuil de variance a ete fixe a 99%. La dimension du voisinage du 
point courant pour le calcul de la coherence est de 25 traces, de 4ms chacune. La zone 
5 reservoir correspond a une epaisseur de 20ms. 

La figure 8 montre les trois campagnes sismiques, ainsi que le cube de coherence 
associe. Les zones montrant les valeurs de coherence les plus faibles semblent se situer a la 
base du reservoir, dans les deux tiers sud. La coincidence entre la situation des puits 
permettant la production et les faibles valeurs de coherence conforte T interpretation en 
10 termes de variations 4D et non pas en termes simplement de bruits sismiques, comme cela 
pourrait etre fait compte tenu de la non parfaite reproductibilite du signal demontree 
precedemment avec V attribut de coherence sur la zone hors reservoir. 

Ceci est confirme avec la figure 9 rassemblant les onze plans temporels du fcube de 
coherence. Meme s'il n'est pas totalement immuable, le tiers nord semble ne pas changer 

15 d'une campagne sismique a l'autre, avec des valeurs de coherence superieures a{ : -0,8 en 
majorite, sur toute Tepaisseur du reservoir : il est a noter tout de meme de legeres 
variations entre les CDP 80 a 90 et les lignes 14 a 20 pour les plans situes 12ms a 16 ms 
sous le toit du reservoir. Le coin sud-est du reservoir reste aussi inchange d'une campagne 
a F autre. Ces zones n'apparaissent done pas etre trop influencees par 1'exploitation du 

20 champ : elles peuvent etre considerees comme une zone reservoir de moins bonne qualite 
en terme de porosite/permeabilite. 

Les larges zones de ties faibles valeurs de coherence a la base coincident avec la 
presence de trois des quatre puits d'injection de vapeur et de recuperation de 1'huile, ainsi 
qu'a la partie sud en dessous de ces puits, ce qui laisserait penser a une invasion de la 
25 vapeur injectee sur cette zone. De meme, la zone de tres faible coherence au toit se situe a 
Taplomb de l'extremite des quatre puits : la aussi, cette zone peut correspondre a une 
remontee de vapeur en bout de puits. 

En revanche, le puits situe le plus au nord coincide avec une zone legerement plus 
coherente au-dela de la ligne 80. Ce puits est situe a la limite avec la zone consideree 
30 comme moins bon reservoir; la vapeur injectee pourrait davantage influencer la partie 
situee plus au sud de ce puits. 
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La figure 10 montre le plan tempore] situe 12ms sous le toit du reservoir extrait des 
cubes de coherence calcules selon les deux autres approches : pour ]a premiere methode en 
prenant en compte une seule variable synthetique (a), deux variables synthetiques (b), pour 
la seconde methode en fixant un seuil de variance a 90% (c), a 95% (d), et a 99% (e), Les 
5 deux cartes (a) et (b) sont tres similaires entre elles, mais egalement tres similaires 
respectivement aux cartes (d) et (e). L'ajout d'une seconde variable synthetique, comme 
pour le cas hors reservoir, ne semble pas changer T interpretation qui pourrait etre faite. 
Globalement, il semble que deux variables synthetiques suffisent a expliquer presque la 
totalite de la variance initiale de chacun des groupes. d'attributs analyses. De meme, la 

10 prise en compte de variance supplemental entre les cartes (d) et (e) ne change pas les 
coherences obtenues, a I'exception de petits details. En revanche, la carte (c) apparait 
beaucoup plus coherente que les deux autres cartes. La part de variance locale 
supplemental prise en compte correspond done, dans la plupart des cas, a une 
information qui est moins commune aux trois cubes consideres. Les valeurs de coherence 

15 obtenues dans ce cas sont plus fortes que les valeurs de coherence obtenues en effectuant 
une moyenne simple (voir la carte correspondante sur la figure 9). 

La figure 11 montre une vue 3D du cube de coherence obtenu avec deux variables 
synthetiques et regroupant les valeurs de coherence strictement inferieures a 0,8. II apparait 
nettement que le tiers au nord est inchange, ainsi que le coin nord-est. De meme, seuls les 
20 deux tiers au sud semblent changer. 

2 - Application a des donnees sismiques avant sommation 

La methodologie peut egalement s'appliquer a la sismique avant sommation : dans 
ce cas, il s'agit de determiner s'il existe des zones coherentes dans les donnees AVO a 
partir de plusieurs cubes sismiques 3D iso-angles ou iso offsets, 

25 Les donnees utilisees sont constituees des cinq cubes iso-angles couvrant un 

reservoir petrolier (chenal avec des depots greseux). L'epaisseur de la sequence etudiee est 
de 38ms. 

La taille de chaque voisinage est de 5 lignes par 5 CDP, soit un total de 25 traces. 
La dimension verticale retenue est de 4 ms, soit trois echantillons temps. Le cube de 
30 coherence a ete calcule selon la premiere methode (moyenne simple) avec un seuil de 
variance fixe a 99%. 
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La figure 12 montre trois des cinq cubes iso-angles utilises (cubes 0°-6° s 12°-18° et 
24°-30°) ainsi que le cube de coherence obtenu. Sur ce dernier cube, les zones les plus 
coherentes apparaissent en orange et en rouge, et les zones les moins coherentes en vert et 
en bleu. Les bordures d'une structure chenalisante y apparaissent nettement sous la forme 
5 de zones coherentes. 

Globalement, les zones les moins coherentes se situent essentiellement dans la 
partie haute de la fenetre reservoir etudiee (figure 13, carte a), a 1'exception d'une petite 
zone tres coherente au nord-ouest correspondant a une forte anomalie d' amplitude se 
retrouvant sur tous les cubes angles. 

10 Dans la partie mediane (carte b), les zones les plus coherentes suivent les contours 

de la forme chenalisante, le chenal en lui-meme correspondant a des valeurs de coherence 
inferieures a 0,8. Dans la partie basse de la fenetre (carte c), les zones incoherentes sont 
moins nombreuses et sont essentiellement situees au nord-est et au sud-ouest. 

Les zones les moins coherentes semblent souligner des zones sismiquement plus 
15 sourdes ou des zones sismiques pour lesquelles les marqueurs ne se retrouvent pas d'un 
cube angle a V autre. 

La figure 14 montre la ligne passant par un puits W2, extraite des cubes sismiques 
0°-6°, 12°-18° et 24°-30°, ainsi que la meme ligne extraite du cube de coherence. Les 
zones correspondant aux chenaux sont relativement bien marquees par des faibles valeurs 
20 de coherence dans leur partie haute", et des valeurs plus elevees dans leur partie basse. Les 
. zones coherentes correspondent bien a des marqueurs de forte amplitude se retrouvant sur 
les differents cubes angles. 

II a egalement ete verifie qu'en diminuant le nombre de traces pris en compte dans 
le voisinage, le cube de coherence obtenu prenait un aspect plus bruite. De meme, il a ete 
25 verifie qu'en augmentant la dimension verticale des traces sismiques du voisinage, le cube 
de coherence obtenu etait lisse verticalement. 

L'attribut de coherence AVO montre done le degre de coherence des cubes 
sismiques extraits dans le voisinage des points et consideres en fonction de Tangle. En 
consequence, les zones incoherentes peuvent etre interpretees soit comme du bruit 
30 sismique, soit comme des fades lithologiques particuliers transparents du point de vue 
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sismique (c'est-a-dire ne montrant pas de reflecteurs), ou comme de fortes variations de 
l'amplitude en fonction de l'angle (dues, par exemple, au contenu fiuide). II est done 
interessant de prendre en corapte cet attribut de coherence dans l'lnterpretation des 
reservoirs, en complement d'autres attributs. 
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REVENDICATIONS 

1) Methode pour mesurer les similarites locales entre un certain nombre p de cubes 
de traces sismiques obtenus par exploration sismique d'un meme volume d'une zone 
souterraine, caracterisee en ce qu'elle comporte les etapes suivantes : 

a) on extrait, de chaque cube de traces sismiques, un voisinage volumique, centre 
sur un meme point courant et compose d'un ensemble de traces sismiques en nombre 
limite ; 

b) on applique la technique d'analyse ACPG a des groupes d'attributs sismiques 
extraits des traces sismiques du dit voisinage volumique pour former des variables 
synthetiques ; 

c) on determine une valeur de coherence a partir des variables synthetiques, que 
Ton affecte au point courant ; 

d) on reitere les etapes a) a c) pour chacun des points communs aux differents 
cubes; et 

e) on regroupe 1 'ensemble des valeurs de coherence pour former un cube de 
coherence traduisant les dites similarites. 

2) Methode selon la revendication 1, caracterisee en ce que, pour chaque point, on 
prend comme valeur de coherence, la valeur moyenne des carres des correlations entre un 
nombre limite k de variables synthetiques et leurs projections sur les cubes dans le 
voisinage du point courant. 

3) Methode selon la revendication 2, caracterisee en ce que Ton determine la valeur 
du nombre k de variables synthetiques comme le plus petit nombre de variables 
synthetiques permettant d'atteindre un seuil de variance explique par les projections des 
variables synthetiques sur chacun des cubes, ce seuil etant prealablement choisi. 

4) Methode selon Tune des revendications precedentes, caracterisee en ce que Ton 
choisit un nombre k de variables synthetiques selon leurs correlations avec les groupes 
d'attributs associes au voisinage volumique du point courant, la valeur de coherence 
attribuee au point courant etant egale a la somme ponderee des carres des correlations entre 
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les variables synthetiques considerees et leurs projections sur les cubes dans le voisinage 
du point courant. 

5) Methode selon la revendication 4, caracterisee en ce que, pour une valeur de 
correlation definie, on choisit comme valeur de ponderation, le pourcentage de variance 

5- explique par la projection de la variable synthetique sur le groupe correspondant divise par 
la somme des variances de toutes les projections des variables synthetiques considerees sur 
le meme groupe. 

6) Methode selon Tune des revendications precedentes, caracterisee en ce que Ton 
fixe un seuil sur le pourcentage de variance expliquee par les projections des variables 

10 synthetiques sur les cubes dans le voisinage du point courant qu'il faut prendre en compte, 
la valeur de coherence etant alors dgale a la somme ponderee des carres des correlations 
entre les variables synthetiques et leurs projections sur les cubes dans le voisinage du point 
courant, de sorte que le nombre de variables synthetiques pris en compte permette 
d'atteindre ce seuil. 

15 7) Methode selon la revendication 6, caracterisee en ce que pour une valeur de 

correlation, on choisit comme valeur de ponderation, p fois le seuil de variance fixe, p etant 
le nombre de cubes sismiques consideres. 

8) Methode selon Tune des revendications precedentes, caracterisee en ce que Ton 
extrait un voisinage volumique de cubes de traces sismiques obtenus a Tissue d'une 

20 campagne sismique 3D, chacun correspondant a un meme angle d'incidence. 

9) Methode selon Tune des revendications 1 a 7, caracterisee en ce que Ton extrait 
un voisinage volumique de cubes de traces sismiques obtenus a Tissue d'une campagne 
sismique 3D, chacun correspondant a un meme deport ou offset. 

10) Methode selon Tune des revendications 1 a 7, caracterisee en ce que Ton extrait 
25 un voisinage volumique de cubes de traces sismiques obtenus par des explorations 

sismiques successives de la zone. 

11) Methode selon Tune des revendications 1 a 7, caracterisee en ce que Ton extrait 
un voisinage volumique de cubes de residus obtenus a Tissu d'une inversion 
stratigraphique avant sommation. 
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12) Methode selon Tune des revendications 1 h 7, caracterisee en ce que Ton extrait 
un voisinage volumique de cube de residus obtenus a l'issu d'une inversion strati graphique 
apres sommation. 

13) Methode selon Tune des revendications 1 a 7, caracterisee en ce que Ton extrait 
un voisinage volumique des cubes de traces inverses avant ou apres sommation, et des 
cubes de residus. 
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